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ABSTRACT 

We present deep submillimeter imaging of two spectroscopically confirmed z ~ 6.5 Lyman-a emitters (LAEs) 
in the Subaru Deep Field. Although we reach the nominal confusion limit at 850/im, neither LAE is detected 
at 850/im nor450/im, thus we conclude that the LAEs do not contain large dust masses (< 2.3 x 10 8 M and 
< 5.7 x 10 8 Mq). The limit on their average Lfir/Luv ratios (< 35) is substantially lower than seen for most 
submillimeter selected galaxies at z ~ 3, and is within the range of values exhibited by Lyman-break galaxies. 
We place upper-limits on their individual star formation rates of < 248 MQyr -1 and < 613 M yr _1 , and on the 
cosmic star formation rate density of the z ~ 6.5 LAE population of < 5.0 x 10~ 2 MQyr _1 Mpc~ 3 . In the two 
submm pointings, we also serendipitously detect seven sources at 850/im that we estimate to lie at 1 < z < 5. 

Subject headings: galaxies:formation - submillimeter - galaxies, starburst - galaxies, high-redshift 



1. INTRODUCTION 

Deep narrow-band imaging at optical wavelengths has 
identified a sizable number of star-forming galaxies at 
z ~ 6.5 through searches for strong Lyman-o: (Lyct) emis- 
sion JAiiki et alJ 120031 Ikodaira et alJl2003t iHu et al] l200l 
Rhoa ds et alJl2004t iTaniguch iet al]l2005l) . At this redshift, 
the universe was less than 10 9 years old and these early sys- 
tems provide insight into many important aspects of the for- 
mation of galaxies at very early times. Although the sample 
size of z > 6 galaxies is increasing through the efforts of large 
surveys, most efforts are focused on optical and near-infrared 
observations and follow-up studies at longer wavelengths re- 
main scarce. 

The detection of z ~ 6.5 star-forming galaxies at submm 
or far-infrared wavelengths would have a number of impor- 
tant implications. Firstly, detection at these wavelengths re- 
quires that these young systems contain substantial amounts 
of dust and, at < 10 9 years after the Big Bang, would pro- 
vide current dust production models with strong constraints. 
Although large dust masses have been found in high red- 
shift (z > 4) quasars and radio ga l axies (e.g. | Q mont e^ alJ 
200 11 iPriddev & McMahor] 12001b [Archibald e t alJ 1200 It 



Robson et al. 2004), quasar winds offer an explanation for the 



production of their dust jElvis. Marengo. & Karovskal2004l) : 
however, such a mechanism cannot be invoked for pure star- 
burst galaxies, i.e. the z ~ 6.5 Lymana emitters (LAEs). Be- 
cause the age of the universe at z ~ 6.5 is comparable to 
the evolutionary timescales of low-mass stars, evolving stars 
cannot be the primary source of dust. Rather, these star- 
forming galaxies can only manufactur e dust via grain for- 
mation in Type II superno vae (SNe) llTodini & Ferrarall2001t 
Morgan & Edmunds 2003); even so, a top-heavy Initial Mass 
Function and very high star formation efficiency also must be 
invoked. 

Secondly, quantifying the amount of dust in galaxies at z ~ 
6.5 is important for measuring their individual star formation 
rates, bolometric luminosities, and by extension, the cosmic 
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star formation rate density (SFRD) at z > 6. The SFRD has 
now been measured for large samples of galaxies to z ~ 5 
and thus far does not exhibit a ny clear sign of turning over 
jLillv et alJ 119951 [Madau et alJl!996t IConnollv et al] 1 19971 
Steide i*eUdIT9^lGiayali^ however, this must 

occur somewhere in the relatively short amount of time be- 
tween z ~ 5 and the epoch of reionization. Currently the ma- 
jority of measurements at z > 6 rely on a relatively small num- 
ber of emission line galaxies fo und through narrow-band op- 
tical searches for Lya emiss ion SHu et all2002t iRhoads et alJ 
120041 Taniguchi et al] l2005l) . While these studies measure a 
SFRD that is apparently an order of magnitude less than at z ~ 
5, their accuracy is hindered by small numbers, uncertainties 
in the completeness corrections, and the unknown amount of 
dust extinction at z > 3. 

The unknown correction for dust at z > 6 is potentially 
substantial, as illustrated by recent submm su rveys (e.g., 
Smail. Ivison. & Blainl f l997 | iLillv et alJ fl999t Webb et al. 
2003b; Chapm an et al J 120051) that show the correction due 
to dust is large at z ~ 3: the dust-uncorrected optical/UV 
estimates of the SFR are one to two orders of magnitude 
lower than the infrared estimates, and the total star for- 
mation in rare, infrared-selected objects is comparable to 
that seen in large po pulations of optically-selected galaxies 
( Chap man et alJl2005l) . Thus, if there exists a population of 
heavily dust-obscured star-forming galaxies at z ~ 6.5, in- 
frared observations are crucial for measuring their total in- 
dividual and integrated star formation rates. 

One might assume that because dust is an efficient ab- 
sorber of Lya photons, the existence of strong Lya emis- 
sion in these systems precludes the presence of significant 
amounts of dust. However, this does not empirically ap- 
pear to be the case: thus far, ~ 30% of the spectroscopic 
redshifts for submm-detected galaxies are, in fact, derived 
from strong Lya in emission ( Chapm an et al J 120051) . Only 
a small number of LAEs at z ~ 2 have been observed in 
the submm, and though none of these have been detected 
(Barger, Cowie, & Sanders 1999; Smail et al. 2003) we note 
that many of the larger, extended Lyman-a blo bs (LABs) have 
been detected in the submm and/or infrared llChapman et alJ 
l200lUGeach et alJ2005UColbert et alJ2 006). 

Here we extend for the first time the search for submm 
emission to two LAEs at z = 6.5. We are aided by the neg- 
ative ^-correction at submm wavelengths that results in equal 
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sensitivity for detecting dust-obscurred star formation at z ~ 6 
as at z ~ 1 ■ The two targeted LAEs were discovered through 
narro w band imaging at A r =9196A in the Subaru Deep Field 
(SDF; Kodaira et al. 2003), and have been spectroscopically 
confirmed to lie at z = 6.541 (SDF-1) and 6.578 (SDF-2). We 
present results at 850 and 450/im for both LAEs; while nei- 
ther object is detected, we attempt to place upper limits on 
the properties of the z ~ 6.5 LAE population using our deep 
submm observations. The paper is organized as follows: in 
§2 we describe the observations and data analysis; §3 presents 
the derived properties of the z ~ 6.5 LAEs; and §4 a brief anal- 
ysis of the serendipitous detections. Our conclusions are in 
§5. We assume an J7m=0.3, 51a=0.7 cosmology with H o =70 
km/s/Mpc throughout. 

2. DATA AND ANALYSIS 
2.1. Observations and Data Reduction 

We observed the objects using the Submilli meter (submm) 
Comm on-User Bolometer Array (SCUBA; Holla ncTet al.l 
1999) on the James Clerk Maxwell Telescope (JCMT) over 
eight partial nights in 2003, 2004, 2005. We obtained data at 
850/im and 450/im simultaneously, but due to beam instabil- 
ities and poor sky transmission the 450/im data are of poorer 
quality than the 850/im data. To fill in the under-sampled sky, 
SCUBA was stepped through a regular 64-point jiggle pattern 
and sky flux was removed to first order by employing a 3 -point 
chop with a chop-throw of 30". The chop position angle was 
held constant in right ascension which results in a characteris- 
tic negative-positive-negative beam shape on the final map for 
real point sources. Telescope pointing was checked every 1.5 
hours, and sky opacity was monitored through sky dips every 
1.5-3.0 hours and in real time along the line of sight using the 
JCMT water vapor meter when operational. 

The data were reduced using a combination of general 
SURF ( Jenness & Lightfoot 1999) routines and custom IDL 
programs written by ourselves. We first applied standard flat- 
fielding and extinction corrections. Although the 3 -point chop 
employed by SCUBA removed sky to first order, residual sky 
flux remained; SURF attempts to remove this through a sub- 
traction of a single median or average sky measurement for 
each second of measurement. To improve on this, we sub- 
tracted a sky plane fit to each second of the data time streams 
and all bolometers, with the bolometers weighted by their in- 
dividual noise. In a single second of integration, any structure 
in the maps is only due to sky flux and thus this does not 
add any systematic offset. A comparison of images produced 
with and without this method indicates a decrease in the noise 
of ^30%. Noise spikes were then iteratively removed using a 
simple standard deviation clipping routine (at 3cr using the en- 
tire array) and the data were regridded onto an astrometrically 
calibrated map using standard Gaussian weighting. 

Total integration times (including the off-source chop time) 
were 47ks for SDF-1 and 19ks for SDF-2. The final images 
have central (unsmoothed) noise values of rms=1.2mJy/beam 
(SDF-1) and rms=1.6mJy/beam (SDF-2); thus both pointings 
reach depths close to the nominal confusion limit at 850/im. 
Note that the noise does not scale perfectly as f" 1 / 2 due to 
weather differences and night to night variations in the noise 
properties of the bolometer array. The jiggle-pattern and ro- 
tation of the Nasmyth-located SCUBA on the sky results in 
uneven coverage of the sky, with the center of the image re- 
ceiving more integration time than the outer edges and hence 
the centers of each image, where the LAEs are located, are 



the deepest points on each map. 

2.2. Source Detection 

To improve point source recovery, we have employed 
the so urce detection and extra ction technique of Eal es et aT] 
(2000) and IWebb et alJ (l2003bl) . In short, the map is itera- 
tively cleaned using a beam template generated from the point 
source calibration maps. Sources are detected by convolving 
the cleaned map with the beam template. Such a technique 
takes advantage of the negative-positive-negative pattern that 
real sources exhibit, but that noise spikes do not, and improves 
the separation of blended sources. 

The noise as a function of position in each map was de- 
termined in the following way. We produced fake bolometer 
time-streams, drawn from the real data using a bootstrapping 
technique. These were passed through our analysis pipeline 
to produce 1000 full map realizations and the noise is taken 
to be the rms variation of each pixel in these maps. This tech- 
nique agrees well with Monte Carlo simulations of the noise, 
assuming Gaussian statistics JWebb et al.l2003bH2005l) . 

2.2.1. SDF-1 & SDF-2 

In Fig. [0 we show the two LAE fields at 850/im after 
smoothing with a 14" Gaussian; the LAEs are at the center 
of each image. Overlaid are the signal-to-noise ratio (S/N) 
contours resulting from the iterative cleaning procedure for 
the 850/im and 450/im maps; the 450/im contours are shown 
over the 850/im grey-scale map to aid in positional compar- 
ison. Neither image contains significant flux at the locations 
of the LAEs. Fig. |2] shows the two 850/im maps after the 
>3.5cr sources have been cleaned and removed, with the LAE 
position marked and the flux at this position listed in Tabled 

Although the maps have a mean flux of zero, there is resid- 
ual positive and negative structure which is not only due to 
noise but also real structure in the extragalactic submm sky 
below the confusion limit of SCUBA. This combination of 
noise and source confusion results in an offset between the 
true and recovered positions of objects; simulations show that 
sources at ~3.5-4u are rarely recovered at their original lo- 
cation and can be shifted up to 8" aw ay (~ 90-95 percentiles 
IWebb et al.l[2003bl llvison et alJl2005l) . Thus, if either of the 
LAEs exhibited significant submm flux, one would expect it 
to be offset from the known position by a few arcseconds; 
however, there are no significant sources within any reason- 
able search radius of the LAEs. 

To determine robust upper limits on the observed flux of the 
LAEs, we placed and recovered sources at the LAE position 
in each map, following the same cleaning procedure described 
above. Fake sources were deemed recovered if they were de- 
tected at >3cr in the cleaned map within 8" of the original 
input position. For the SDF-1 field, this analysis indicated a 
detection threshold that was in good agreement with the limit 
expected from the noise map. Although the LAE position is 
coincident with the off-beam of a bright (3.8mJy) source, we 
were able to separate the contributions from this and the input 
(fake) source and recover the fake source at > 3er within 8" 
of the input position, provided the input flux was >1.5mJy. 

The SDF-2 map is shallower than SDF-1 and appears to 
suffer from greater source confusion in its central region. Al- 
though the LAE is not confused with a bright neighbor, struc- 
ture beneath the 3er detection limit makes it difficult to recover 
flux at the LAE position. Specifically, it lies in a negative re- 
gion of the map (which is not statistically significant), and 
although this minimum could be the off-beam of one or more 
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faint neighbors, they lie below the confusion limit and thus 
cannot be cleaned from the map. Our simulations indicate 
that we are unable to recover sources at this position above 
3cr below an input flux of 3.8mJy to within 8" of the input po- 
sition. These recovery limits are listed in Table[0for 850/im 
and 450/im where both were determined in the manner de- 
scribed above. 

2.2.2. Serendipitous Detections 

Although the z ~ 6.5 LAEs are not detected at submm 
wavelengths, the deep pointings reveal multiple sources near 
the LAEs; the sources detected in the two fields are listed in 
Table |3 SDF-1 contains five sources detected above 3.5ct 
at 850/im while SDF-2 contains two such sources. Of these 
seven 850/im sources, two are also cleanly detected at 450/im 
with S/N> 3.5 and positional offsets of 4.0". A third 850/rni 
source has two possible 450/im counterparts: a 4. Oct detec- 
tion 10" from the 850/im position, or a 2.8ct peak offset by 
5.3". Because a 10" offset between the 450/im and 850/im 
positions is statistically unlikely, we assume that the closer 
2.8ct peak corresponds to the 450/im emission for this source; 
however, we note that none of our final conclusions are greatly 
altered if the more distant object was the correct counterpart. 
Also listed in Table |2] are three additional sources detected 
only at 450/rni (all at >3.5ct). 

3. LYMAN-a EMITTING GALAXIES AT Z ~ 6.5 

3.1. Star Formation Rates 

A galaxy's far-infrared luminosity is due to thermally ra- 
diating dust and is correlated with the galaxy's current rate 
of star formation, provided the energy contributi on from an 
active gala ctic nucleus (AGN) is negligible (Kennicutt 1998; 
iBelll 12003). Locally, the most intense starburst galaxies are 
also the most luminous in the infrared; however, these systems 
do not show a simple relati on between infrared ex cess and UV 
slope (the IRX-/3 relation; Goldader et al. 2002). Therefore 
extrapolating from the UV to a total Lboi and corresponding 
star formation rate is not trivial, and in many cases it may not 
even be possible. 

This situation is mirrored at high redshift (z ~ 2-3) where 
luminous (L > 10 u L m ), dusty starbursts are m ore preva- 
lent (e.g., iBorvs et all 120021 [Webb et al. 2003b). For such 
systems, the rest-frame UV estimated SFR can be two or- 
ders of magnitude less than that im plied by the rest-frame far- 
infrared emission (Chap man et al.ll2~005l) . Far-infrared imag- 
ing offers an independent and, for the most luminous systems, 
a more accurate measure of their total star formation rates and 
bolometric luminosities. Moreover, because of the flat flux- 
redshift relation for infrared luminous galaxies at 850/im be- 
yond z ~ 0.5, such systems are equally detectable at z ~ 6 
as they are at z ~ 1, unlike optical and near-infrared observa- 
tions. 

However, converting infrared luminosities to star formation 
rates suffers from a number of uncertainties. While the em- 
pirically measured local SFR-Lir relation does hold over four 
orders of magnitude in Lir, the scatter about this relation is 
approximately a factor of two. At high redshifts, the scatter 
could be much greater and/or systematically offset from the 
local relation as unce rtainties due to diff erences in the Initial 
Mass Function (IMF; Baugh et al. 2005), differences in dust 
properties such as temperature or emissivity, and the contam- 
ination from AGN increase. This is especially relevant at the 
redshifts considered here where we are beginning to probe the 
very first generation of stars and dust. 



With these caveats in mind, we have employed the follow- 
ing method to convert our limit on the observed Ssso/im to a 
star formation rate. We adopt a modified black-body spectral 
energy distribution (SED) with a temperature of 36K, a dust 
emissi vity of /?= 1 .5. Following Blain. Barn ard. & Cha pman 
(2003) we apply a Wien correction short-ward of 80/im with 
a slope of a=2.9. These parameters are chosen to match 
the measureme nts of the sub millimeter luminous popula- 
tion at z ~ 3 (Chapman et all I2005I) . but they are also in 
good a greement wi th the results from local ULIRG ga laxies 
[Dunne et al]l2000t iBlain. Barnard & Chapmanll2003l) . Us- 
ing this SED, we calculate Lir (i.e., L(8- 1000/im)) for the 
two LAEs and use the relation from Bell (2003) to convert 
this to a star formation rate. The results are included in Ta- 
ble|2] Differences in the assumed temperature of ^30% (the 
inter-quartile range of Chapman et al.) result in a difference in 
the inferred star formation rates of ~70%, in the sense that as- 
suming hotter dust results in larger a infrared luminosity and 
therefore a larger inferred star formation rate. 

These SFR limits can be directly compared to those de- 
termined through optical measurements. The average dust- 
uncorrected Lya-estimated rate for the two galaxies is ^7 
M Q yr -1 , and that estimated from the UV continuum is slightly 
higher ~22 M Q yr _1 (Kodaira et al. 2003). While our average 
upper-limit of ~ 430 MQyr" 1 is an order of magnitude larger 
than this estimate, it is still a useful first constraint on the in- 
frared properties of these systems; recall the results of the 
submm and radio studies at z ~ 3 that reveal star formation 
rates in excess of ^1000 M^yr " 1 in systems for whi ch the 
UV prediction is ~ 10 M yr _1 dChapman et alJl2005l) . The 
two LAEs discussed here are definitely not members of such 
a luminous population. 

3.2. Dust Content 

Our deep submm imaging offers a unique opportunity to 
constrain the mass of dust present in early starburst galaxies, 
a measurement that has thus far only bee n possible at these 
redshifts for infrar ed luminous quasars ( Isaa k et all 2002; 
Pridd evet alj|20 03). Using the same parameters we adopted 
in §3.1 (T = 36K, = 1 .5), we calculate the dust mass follow- 
ing iHildebrandlllSsl : 



M d : 



(\+z)k v B v (J) 



(1) 



where k v is the dust emissivity, v is the rest-frame frequency, 
and B, y is the Plan ck function at a temperature T. Following 
iHildebrandl 11983). we adopt «;(125/im) = 1.875 m kg and 
scale with wavelength by n v oc v$ . The results (see Table|2j 
are consistent with the dust mass of the lensed Lyman-break 
galaxy MS 1512-cB58 which is roughly half the value of our 
upper limits for similar dust parameters ( Sawicki 2001). 

Thus far, quasars are the only submm-detected systems that 
are confirmed to lie at z ~ 5-6 and SED arguments suggest 
that their submm luminosity is dominated by thermally emit- 
ting dust rather than by non-thermal emission from the AGN 
dPriddev et alj |2003). The dust mass estimates of the quasars 
range from a factor of 2 to 10 higher than our 3ct upper limits 
(for the same dust parameters). Unless the dust temperature 
of the LAEs is much hotter than for AGN, we would have de- 
tected objects with comparable dust masses in our data. The 
lower dust content of the tw o LAEs is consistent with the re- 
cent mid-infrared results of Char y. Stern. Eisenhardtl (2005) 
who report evidence for A v = 1.0 mag of dust extinction in 
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a single z ~ 6.5 LAE, a v alue that is in the range of normal 
galaxies (Choi et al. 200a). 

While a positive submm detection for either LAE would 
have provided dust formation models with further empiri- 
cal constraints, the non-detections discussed here are unfor- 
tunately ambiguous. The existence of large dust masses in 
these LAEs would raise the important question of how large 
amounts of dust can form in the time between the formation 
of the first objects and z ~ 6.5. Although the two z ~ 6.5 
LAEs have distinctly less dust than the submm detected z > 5 
quasars, the z ~ 2-3 SCUBA population, and many of the z ~ 
2-3 Lya blobs, we cannot constrain them to be less dusty than 
present-day normal starburst galaxies. 

The lack of a detection for either systems raises the fol- 
lowing question: are the z ~ 6.5 LAEs even old enough 
to have built up a sizeable dust reservoir through stellar 
evolution? Perhaps large quantities of dust exists only 
in the oldest systems at z > 6 and/or those for which 
non-stellar dust production mechanisms are available: i.e., 
quasars. The Lya equivalent-widths of the LAEs of 100- 
300A (with uncertainties of > 50%) imply ages of < 10 7 
years (Ma lhotra & Rhoadsll2002t ICharlot & Fal]lll993lh and 
this is too young an age for large amounts of dust to have 
been produced by any known mechanism. Low mass stars 
have not had sufficient time to evolve and even supernovae 
models have difficulty p roducing dust in such young objects 
(Todi ni & Ferraral200ll) . 

Still, such large Lya equivalent widths have been measured 
for members of the submm detected population at z ~ 2-3 
(Chap man et al!2003l) . These systems clearly have large dust 
masses but their Lya EWs imply similar ages to the LAEs. 
Thus while the Lya EW may be a reasonable indicator of the 
age of the current starburst, it cannot provide complete infor- 
mation on the total star formation and dust production history 
of a galaxy. If the star formation histories of the z ~ 6.5 LAEs 
are extended over the first Gyr of the universe significant dust 
mass cannot be ruled out on the basis of the age of the current 
starburst. 

3.3. Comparison to Lower Redshift Populations 

Using the submm upper-limit on the LAEs of Lfir < 3 x 
10 12 L (T) results in an upp er- limit on the Lfir/Luv ratio of 
< 35 (iKodairaet alH 2003). Our limit on the Lfir/Luv ra- 
tio is substantially lower than the values of ~ 1000 that 
are observed for the submm-luminous population at z ~ 3 
(Chap man et al.ll2~005l) . Rather, a Lfir/Luv ratio of < 35 is 
within the range of v alues exhibited by Lyman-break galax- 
ies of Lbtp/Lttw < 10 dAdelberger & Steidell2000tlWebb et all 
I2003al Goldade reral.1120021) ; not e that at z ~ 3, LAEs make 
up ~ 30% of the LBG population dShaplev et alJ2003l) . 

In contrast to the LAEs, many of the extended Lya emit- 
ters at lower redshifts, the so-called Lya Blobs (LABs), are 
infrared bright and have been detected by SCUBA or, more 
recently, at 24/im by the Spitzer space telescope. In Fig. [3] we 
compare the Lya luminosities to the bolometeric luminosities 
(assuming Lir ~ Lboi) of the two z ~ 6.5 LAEs and the lower 
redshi ft, in frared detected LAEs , as discussed in lGeach et al.l 
(2005) and IColbert et alJ (120061) . We divide the points into 
highly extended (diameter> 50 kpc) and compact (diameter< 
50 kpc) objects. There may be a weak tendency for the most 
extended and Lya luminous LABs to have the highest infrared 
inferred bolometric luminosities, though note that two of the 
LABs with LLy Q < 10 43 erg s" 1 have Lboi comparable to the 
most extended and luminous LABs. The z ~ 6.5 LAEs dis- 



cussed here are small (< 4 kpc) with LL yQ comparable to the 
smaller IR bright LABs, but their Lboi are at least an order of 
magnitude lower. 

A difference in the Lboi of different Lya emitting systems, 
assumes however that all of the objects have the same dust 
temperature, and given the large differences in the extent of 
their Lya emission this may not be the case. If the spatial ex- 
tent of the dust scales with that of the Lya emission the tem- 
perature of the dust may also vary, such that the far-infrared 
flux of larger Lya systems is dominated by cold diffuse dust, 
while the compact systems are dominated by similarly com- 
pact hot dust. Given the same intrinsic source luminosity, 
850/im measurements are biased toward cold dust and thus 
would not detect the hotter compact systems; that is, adopting 
a temperature which scales with source size could reduce or 
remove the trend in Fig|3]for the 850/im measurements. How- 
ever, while this might explain the trend seen for the 850/im 
detected objects, such an effect cannot account for the prop- 
erties of the 24/im detected systems, which appear to follow 
the same trend. For these objects Lboi is calculated from the 
rest-frame mid-infrared which does not suffer from the same 
bias as the 850/im data, and thus the smaller LAEs may truly 
be lower luminosity systems. 

An alternative explanation for the range in Lboi lies in the 
physics driving the intense energy production. Locally, lumi- 
nous and ultraluminous infrared galaxies (U/LIRGs) are set 
apart from normal galaxies by their high merger rate and their 
ten dency to host AGN, especi ally at the high luminosity end 
(see Sanders & Mirabel ( 1996) for a review). In these systems 
a major gas rich merger appears to be th e trigger for both 
star formation and AGN activ ity (e.g. iMihos & Her nauist 
1996; Barnes & Hernauist 1996). The preponderance of AGN 
and/or merger signatures persist in infrare d luminous sys- 
tems at high-redshift such as the SMG s ( Alexand er et alJ 
l2003l [Webb et aTl l2003bt IConselice et afll2003l) and some of 
the galaxies as sociat ed with LABs (e.g., |Sm ail et al.l y222f 
Chap man etaLll2005t IWeidinger etalll2005t iDev et alJl2005t 
IColbert et alJ200 6). Infrared faint systems on the other hand, 
such as the SDF LAEs and the bulk o f the LBG population 
(IWebb et alJl2003at iHuang et aUl2005l) do not as a popula- 
tion, exhibit such propert ies dSteidel et alJl!996tlHuang et al.l 
120051 tLehmer et al. 2005). This has lead many authors to con- 
clude that the violent infrared luminous phase of galaxy evo- 
lution is induced through gas rich mergers, and thus systems 
such as the LAEs studied here are infrared-faint because they 
are not currently experiencing such an event. 

3.4. Star Formation Rate Density 

Using our SFR upper-limit and the luminosity fun ction of 
LAEs at z ~ 6.5 presented by Kashikaw a et al.l (120061) . we can 
estimate the maximum contribution by these objects to the to- 
tal SFRD at this redshift. The two LAEs discussed here have 
< L(Lya) > = 7.5 x 10 42 ergs/s and < SFR >< 430 Mgyr" 1 . 
At this luminosity, they represent the brighter end of the l umi- 
nosity function as published by Kashik awa et al.l (120061) and 
have a cumulative number density of ~ 2xl0~ 5 Mpc" 3 . If 
we take the SFR upper-limit determined here to be represen- 
tative of the entire LAE popu lation above L (Lya) = 2x 10 42 
ergs/s, the depth of the Taniguchi et al. (2005) sample, we cal- 
culate their contribution to the SFRD of to be < 0.05 M Q yr _1 
Mpc" 3 . We remind the reader that this limit is based on a 
3cr upper-limit on the flux, and could systematically shift by 
^70% for a dust temperature range of ±30%. If we assume 
instead that the SFR scales directly with L(Lya) and integrate 
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over the Lya luminosity function of Kashik awa et al.l J2006I) 
to L(Lya) = 2x 10 42 ergs/s, this changes the upper-limit to < 
0.03-0.3M Q yr _1 Mpc , where the range corresponds to the 
range of Schechter parameters describing the number density. 

In Fig. |4] we compare our submm limits on the SFRD to 
measurements from other LAE studies at z > 3 as well as the 
submm population at z ~ 2 - 3 (see figure caption for refer- 
ences). The new z ~ 6.5 limits are comparable to that mea- 
sured for the dust-obscurred population at z ~ 2 - 3 and thus 
does not constrain any turnover. However, we stress that these 
are upper limits on the SFRD from only two LAEs at z ~ 6.5. 
Expanding our submm survey to encompass the current sam- 
ple of z ~ 6.5 LAE candidates in the SDF (58) and, e.g. stack- 
ing the objects, would increase our depth by a factor of 5, and 
therefore begin to confirm a decline in the infrared-measured 
SFRD from z ~ 3 to z ~ 6. 

4. SERENDIPITOUS SUBMM SOURCES 

Although no detection is made of the LAEs themselves, 
the submm maps do contain a number of other sources (see 
Table |2j. SDF-1 in particular contains five objects detected 
above 3.5er at 850/rni, all of which are within 1 arcmin of 
the LAE. This represents a modest excess over the number 
counts expected from confusion limited blank field surveys of 
z < 2 sources (e.g., Webb et al. 2003b). Whether the submm 
sources are at the same redshift or possibly even at z ~ 6.5 are 
intriguing possibilities. 

At lower redshifts of z ~ 3, there is tentative evidence from 
pair counts and correlation analyses that the submm-luminous 
population is strongly clustered; indeed, targeted imaging of 
rare, high-redshift systems such as high-redshift radio galax- 
ies (Qtevens et al. 2003) have revealed over-densities of sub- 
millimeter sources. To explore the possibility of whether the 
submm sources near the LAEs belong to a single massive 
structure, either at the redshift of the two LAEs or in the fore- 
ground, we attempt to constrain their redshifts. 

Two of the five 850/im sources in the SDF-1 pointing and 
one in the SDF-2 pointing are also detected at 450/im, and 
their submm flux ratios can be used as a very rough indication 
of their redshift (Fig. [5}- Based on these ratios, it is unlikely 
that the 450/im-detected objects lie at the same redshift as the 
z ~ 6.5 LAEs, unless they have extremely hot dust (^80K). 
Given a reasonable temperature range, the submm sources lie 
between l< z <5. The upper-limits on the four objects not 
detected at 450/im are consistent with this lower limit of z > 
1. 

Considering the large uncertainties on the 450/im-850/im 
flux ratios, it is possible that all of the submm sources lie 
at the same foreground redshift, however, we cannot con- 
firm nor deny this with the current data-set. We have at- 
tempted to identify the optical counterparts to these objects 
using the deep optical imaging available on the SDF, but the 
number density of optical galaxies is simply too great to al- 
low unambiguous iden tifications. Such a n exercise requires 
deep Spitzer imaging (fAshbv et alJ 120061) or radio observa- 
tions Jlvison et a l. 2002). We note that Spitzer imaging of the 
SDF will become available in the next year, at which point we 
will attempt to identify the optical/near-IR/mid-IR counter- 
parts to the submm sources to better constrain their redshifts. 

5. CONCLUSIONS 

Quantifying the dust content in starburst galaxies at z > 6 is 
important both for determining cosmological properties such 
as the star formation rate density and for testing stellar evo- 



lution models, in particular how quickly large dust masses 
can form. Using deep submm observations of two spectro- 
scopically confirmed Lyman-a emitting galaxies at z ~ 6.5 in 
the SDF, we show that these starburst galaxies do not contain 
large quantities of dust (< 4x 10 8 M Q ). Our Sgso/im imaging 
reaches the confusion limit and enables us to place upper lim- 
its on their star formation rates of < 430 M yr" 1 . The two 
LAEs are markedly less infrared luminous and less dusty than 
the submm detected population at z ~ 2-3, the extended Lya 
Blobs at z ~ 2-3, and the z > 5 submm detected quasars. 
Rather, the LAEs have Lir/Luv ratios of < 35 that are more 
in line with the Lyman-break population at lower redshifts. 

The submm observations at z < 3 provide a measurement of 
the star formation rate density at z ~ 6.5 that is independent of 
the UV determined value. Assuming the two LAEs are rep- 
resentative of the entire z ~ 6.5 L AE population and using 
the Lya luminosity function from Kashikawa et al. (2006), 
we can place an upper limit on the LAE contribution to the 
dust enshrouded SFRD at z = 6.5 of < 0.05 M Q yr" 1 Mpc" 3 . 
Although consistent with submm measures of the SFRD at 
lower redshifts, and therefore unable to constrain a turnover, 
we stress that this value is extrapolated from only two objects. 
More observations are clearly needed. 

We recognize that because our results are based on only two 
z ~ 6.5 LAEs, they are tentative. However, our observations 
have yielded the first constraints on the dust content and total 
SF rates of such systems, and enable us to begin comparing 
the z ~ 6.5 LAE galaxies to the growing number of z > 5 ob- 
jects to better understand how these early populations overlap. 
Our deep submm imaging has also revealed a large number of 
serendipitously detected faint submm sources that are likely to 
be at 1 < z < 5. A follow-up study on these objects can pro- 
vide much needed information on the faint (Ssso^m < 8mJy) 
submm population. 

Further constraints on dust enshrouded star formation at 
very high redshifts are necessary and can be obtained through 
infrared observations of larger samples of high-redshift LAEs. 
Indeed, our program was originally intended to target a larger 
number of z ~ 6.5 LAEs, but it was cut short due to SCUBA's 
unfortunate demise. With the commissioning of SCUBA-2, 
such a study will again become possible and will provide 
a more statistically representative sample. In parallel, con- 
straining the counterparts of submillimeter-selected galaxies 
that do not have radio counterparts in deep radio data will 
quantify the extent of the high-redshift {z > 3) tail of the sub- 
millimeter population. Direct detections of dusty star forma- 
tion at high-redshift, to the level of a few lOMQyr" 1 (as are 
currently predicted in the optical), also will be possible in the 
future with the Atacama Large Millimeter Array. 
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Table 1. Submillimeter measurements" of Lyq-emitters 



Name 




S850 M m(mJy) 


S85u,im(mJy) 


S450fim (mJy) 


SFR IR (M yr- 1 ) b 


SFRL^OVlgyr" 1 ) 


M d 1O 8 M 






measured 


limit 


limit 








SDFJ132415.7+273058 (SDF-1) 


6.541 


-0.027 


< 1.5 


< 9.0 


< 248 


9.1 ±0.8 


< 2.3 


SDFJ132418.3+271455 (SDF-2) 


6.578 


-1.0 


< 3.8 


< 9.5 


< 613 


5.1 ±0.2 


< 5.7 



a 3<r upper limits are given; see §2.2.1 for discussion of the method. 
b See §3.1 for description of the st ar form ation rate estimate. 
c Values taken from Ko daira et alH2003T) 



Table 2. Source Catalogs for the two LAE fields 



Name RA (J2000) Dec (J2000) S 8 so M m (mJy) S/N-850/zm S 4 5 0Aim (mJy) S/N-450A*m A position 

850/im-450Aim (") 



SDF1-850-1" 


13:24:17.80 


27:30:41.5 


5.1 


7.4 


11.9 


2.8 


5.3 


SDF1-850-2 


13:24:14.08 


27:31:06.5 


3.8 


6.5 


12.8 


3.5 


4.0 


SDF1-850-3 


13:24:17.77 


27:31:33.5 


4.6 


6.3 


< 11.4 






SDF1-850-4 


13:24:15.21 


27:31:34.0 


3.1 


4.8 


< 12.0 






SDF1-850-5 


13:24:16.38 


27:30:00.0 


3.0 


4.0 


< 12.6 






SDF 1-450-1 


13:24:18.33 


27:30:34.5 


< 2.0 




15.1 


4.1 


10.0 


SDF 1-450-2 


13:24:18.25 


27:30:16.0 


< 2.2 




16.8 


4.3 




SDF2-850-1 


13:24:20.46 


27:15:28.8 


5.5 


5.5 


18.8 


4.4 


4.0 


SDF2-850-2 


13:24:13.44 


27:15:09.8 


4.1 


3.9 


< 30.3 






SDF2-450-1 


13:24:15.63 


27:14:42.2 


<2.8 




16.0 


3.5 





"There are two possible 450/im counterparts to this source. The first, which is listed here, is a 2.8<r peak at an offset of ~5" from the 850/^m position. The 
second possibility is source SDF1-450-1 (also listed) which is a >4.0cr detection but at a 10" offset. Please see the text §2.2 for discussion. 
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FIG. 1 . — The submm maps of the two LAE fields. The grey scale corresponds to flux density, with white denoting positive flux, smoothed with a 14" Gaussian. 
All four images correspond to 850/im; the top two panels show SDF-1 and the bottom two show SDF-2. Overlaid on the left two panels (a,c) are the 850/im S/N 
contours, determined through our iterative cleaning algorithm. The right two panels show 450/im contours overlaid on the 850/im maps to illustrate positional 
coincidence between the two wavelengths. In all four panels, the contour levels start at 3<r and increase in lcr steps. East is to the left and north is to the top. 
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FIG. 2. — The submm images centered on each LAE, after the iterative cleaning has removed all >3.5<r sources: the left panel corresponds to SDF-1 and the 
right panel to SDF-2. The LAEs lie at (0,0) in both images (denoted by the intersection of the dotted lines), and neither are detected at 850^im or 450/jm. The 
measured flux at (0,0) in each image is listed in TablelTI East is to the left and north is to the top. 
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FIG. 3. — The correlation between the flux in the Lyct line and the bolometric luminosity as inferred from single IR measurements at 24/^m or 850um_for 
extend ed Lya emitters (LABs) at 2 < z < 3, compared to the two SDF LAEs. For consistency we have here adopted the same SED parameters as Ge ach et all 
(2005) of T=40K, a=4.5, /3=1.7; however the parameters are similar to those discussed in the text and do not result in a significant change in L|, i . The dashed 
line is not a fit to the data but corresponds to the case where L(Lya) = 0.001 L^ol 1 Geach et al. 2005). Open points denote highly extended objects with sizes > 
50 kpc 2 , and solid points denote objects that are more compact than this. 
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4 F i i i i i i " 

1 2 3 4 5 6 7 

redshift 

FIG. 4. — The cosmic comoving star formation rate density (SFRD) as a function of redshift. For clarity we show only two populations on this plot which are 
most relevant to this work: the submm luminous galaxies (shaded area) and the LAEs (all other points). The lower bound on the submm shaded area corresponds 
to the measured SFRD from Sg5Q ura > 5mJy galaxies and the upper boun d shows the exp ecte d values when the population is corrected for completeness down to 
S850um = lmJy. Points are Aiiki et al. 1 2003, upside-down triange), Rhoads et al. 1 2003, star) Cowie & Hu 1 1998, triangle), Kudritzki et al. 1 2000, open square), 
and Fuiita et al. 1 2003, diamond). The lower open circle corresponds to the LAE population of Taniauchi et al. 1 2005) from which the two LAEs studied here 
are drawn; and the upper open circle shows the r evised SFRD of this population based on the dust extinction measured by Charv. Stern, Eisenhardl 1 2005). Our 
upper-limit, based on the Taniguchi et al. 1 2005 ) number counts is shown by the large filled circle, and the range in upper-limits calculated from the Schecter- 
function fits of Kashikawa et al. 1 2006) and assuming a linear relation between Lya line strength and IR-estimated SFR is shown by the two smaller solid points. 
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0.1 1.0 
redshift 

FIG. 5. — The 450/xm-850/jm flux ratio as a function of redshift. The three curves correspond to three simple spectral energy distributions: grey body curves of 
varying dust temperature and an single dust emissivity of /3= 1.5. The three horizontal lines show the ratio for the three 850/im and 450^tm detected sources, with 
the single error-bar in the lower left showing the uncertainty in the ratio (approximately equal for all three objects). The four solid stars denote 3<r upper limits 
for the 850/im sources with no 450/xm detection, arbitrarily placed at z = 6 for clarity. Shown for comparison are the measured ratios (diamonds) or upper-limits 
(circl es) for 850/jm-selected sources with spectroscopic redshifts i Scott et al. 2002; Webb et al. 2003c; Clements et al. 2004; Borvs et al. 2004; Chat>man et al. 



